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Rothia  mucilaginosa  is  an  opportunistic  pathogen  in  the  human  oral  cavity  and  pharynx.  We  found  that  R.  mucilaginosa  DY-18, 
a  clinical  isolate  from  a  persistent  apical  periodontitis  lesion,  had  biofilm-like  structures.  Similar  structures  were  also  observed 
on  R.  mucilaginosa  ATCC25296.  To  further  study  these  structures,  we  determined  the  complete  genome  sequence  of  DY-18  and 
found  it  a  2.26-Mb  chromosome.  Regarding  stress  responsive  systems  known  to  affect  biofilm  formation  in  many  bacteria,  DY-18 
genome  possessed  only  two  sigma  factor  genes.  One  of  these  encoded  an  additional  sigma  factor  whose  promoter-binding  activity 
may  be  regulated  in  response  to  environmental  stimuli.  Additionally,  several  genes  assigned  to  two-component  signal  transduction 
systems  were  presented  in  this  genome.  To  the  best  of  our  knowledge,  this  is  the  first  complete  genome  of  R.  mucilaginosa  species 
and  our  data  raise  the  possibility  that  this  organism  regulates  the  biofilm  phenotype  through  these  stress  responsive  systems. 


1.  Introduction 

R.  mucilaginosa,  formerly  classified  as  Stomatococcus  muci- 
laginosus  [1],  is  a  gram-positive,  coagulase-negative,  non- 
spore-forming  coccus  considered  as  a  part  of  the  commensal 
flora  of  the  oral  cavity  and  the  upper  respiratory  tract 
in  humans  [2].  There  have  been  an  increasing  number  of 
infections  reported  for  R.  mucilaginosa,  especially  among 
immunocompromised  patients  [3,  4].  Therefore,  this  organ¬ 
ism  has  received  attention  recently  as  a  potential  emerging 
opportunistic  pathogen.  This  organism  is  also  known  to  act 
as  an  initial  colonizer  on  the  tooth  surface  in  the  oral  cavity 
[5] .  Contrarily,  several  studies  suggested  that  R.  mucilaginosa 
belongs  to  the  healthy  oral  microbial  community  [6]  or 
constitutes  a  microbial  community  whose  productivity  of 
malodor  is  lower  than  those  of  other  communities  without 
R.  mucilaginosa  [7].  So,  the  factors  making  a  divide  for  this 
organism  as  an  opportunistic  pathogen  or  as  a  member 


of  the  healthy  oral  microbiota  are  still  obscure.  A  clinical 
strain,  R.  mucilaginosa  strain  DY-18  (DY-18),  was  isolated 
from  the  bacteria  that  were  present  in  a  persistent  apical 
periodontitis  lesion  [8].  Apical  periodontitis  is  a  relatively 
common  inflammatory  disease  in  dentistry,  and  a  wide 
variety  of  bacterial  genera  have  been  implicated  as  putative 
pathogens  [9-11].  If  allowed  to  progress,  apical  periodontitis 
can  cause  the  destruction  of  supporting  connective  tissues 
and  bones,  ultimately  resulting  in  tooth  loss.  The  ability 
to  produce  extracellular  matrix  and  to  form  biofilms  is 
recently  considered  to  be  crucial  for  microorganisms  that  are 
present  in  a  root  canal  to  resist  the  intraroot  canal  procedures 
of  disinfection,  to  occupy  apical  foramina  of  teeth,  and 
to  cause  persistent  chronic  inflammatory  lesions  [9,  12- 
14].  For  future  studies  to  genetically  clarify  the  mechanism 
of  matrix  production,  we  have  determined  a  complete 
genome  sequence  of  DY-18  with  a  combined  strategy  of 
pyrosequencing  and  paired-end  sequencing. 
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Figure  1:  Scanning  electron  micrographs  of  Rothia  mucilaginosa  strains  DY-18  (a)  and  ATCC25296  (b)  showing  dense  meshwork-like 
structures  around  the  cells  (bar  =  3  ^m). 


2.  Methods 

2.1.  Bacterial  Strains  and  Culture  Conditions.  DY-18  was 
from  our  stock  culture  collection  isolated  from  a  persistent 
apical  periodontitis  lesion,  where  the  presence  of  bacteria 
was  noted  in  root  canals  after  numerous  treatments  [8].  R. 
mucilaginosa  ATCC25296  (ATCC25296),  a  reference  strain 
for  R.  mucilaginosa,  was  obtained  from  American  Type 
Culture  Collection  (Manassas,  VA,  USA).  All  strains  were 
grown  on  blood  agar  plates  (BBL  Microbiology  Systems, 
Cockeysville,  MD,  USA)  at  37°  C  for  24  hours  aerobically. 

2.2.  Scanning  Electron  Microscopy.  To  examine  cell  surface 
structures,  scanning  electron  microscopy  (SEM)  was  per¬ 
formed.  Bacteria  grown  on  blood  agar  plates  for  24  hours 
were  collected  on  a  piece  of  filter  paper  (Glass  fiber  GA55, 
Toyo  Roshi,  Tokyo,  Japan),  fixed  with  2%  glutaraldehyde 
in  0.1  M  phosphate  buffer  for  2  hours  and  1%  OSO4  in 
0.1  M  phosphate  buffer  for  1  hour  at  4°C,  and  dehydrated 
through  an  ethanol  series  and  2-methyl-2-propanol  followed 
by  platinum  ion  coating  (E-1030,  Elitachi,  Tokyo,  Japan). 
Specimens  were  examined  with  a  scanning  electron  micro¬ 
scope  (S-4800,  Hitachi)  at  an  accelerating  voltage  of  3  kV. 

2.3.  Genome  DNA  Extraction.  To  extract  genomic  DNA,  DY- 
18  was  grown  in  a  trypticase  soy  broth  (BBL  microbiol¬ 
ogy  systems)  with  0.5%  yeast  extract  (Difco  Laboratories, 
Detroit,  MI,  USA)  for  24  hours.  Genomic  DNA  was  extracted 
using  MagExtractor  (Toyobo,  Osaka,  Japan)  according  to  the 
manufacture’s  protocol. 

2.4.  Whole  Genome  Sequencing.  The  complete  genome 
sequence  of  DY-18  was  determined  using  a  combination 
of  454  pyrosequencing  and  paired-end  sequencing  [15] 
performed  using  the  Roche  Genome  Sequencer  FLX  Systems 
(Roche  Diagnostics,  Tokyo,  Japan)  and  Applied  Biosystem 
3730  DNA  Analyzer  (Applied  Biosystems,  Foster,  CA,  USA). 
The  sequence  data  were  assembled  with  GS  de  novo  Assem¬ 
bler  Software  (version  1.1.03.24)  (Roche  Diagnostics).  The 
putative  gaps  between  contigs  were  closed  by  primer  walking 


Table  1:  The  allocation  of  gene  function. 


Functional  category 

Number  of  genes 

Transcription 

122 

Translation,  ribosomal  structure  and  biogenesis 

158 

Energy  production  and  conversion 

114 

Carbohydrate  transport  and  metabolism 

117 

Amino  acid  transport  and  metabolism 

207 

Nucleotide  transport  and  metabolism 

61 

Coenzyme  transport  and  metabolism 

100 

Lipid  transport  and  metabolism 

66 

Signal  transduction  mechanisms 

72 

and  sequencing  of  PCR  amplicon  over  the  gaps.  Sanger 
reads  were  incorporated  into  the  assembly  using  Sequencher 
(version  4.7,  Gene  Codes,  Ann  Arbor,  MI,  USA).  Contig 
editing  and  primer  design  were  performed  using  the  in  silico 
Molecular  Cloning  Genomics  Edition  (version  4.1,  in  silico 
biology,  Yokohama,  Japan).  Bases  in  the  assembled  contigs 
carried  a  Phred- equivalent  quality  score  of  40  or  above  (error 
rate  ^  0.010%). 

The  sequence  was  annotated  using  the  coding  region 
identification  tool  invoking  comparative  analysis  (CRITIC A) 
(version  1.05b)  and  the  Gene  Locator  and  Interpolated 
Markov  ModelER  (Glimmer)  2  (version  2.10). 

2.5.  Bioinformatics  Analysis.  Since  sigma  factors  play  a 
critical  role  in  regulating  biofilm  formation  in  some  bacteria 
[16-20],  we  also  attempted  to  identify  homologs  of  sigma 
factors  on  the  genome.  Amino  acid  sequence  of  Escherichia 
coli  K-12  RpoD  (GenBank  ID  (GI):  1789448),  RpoE  (GI: 
1788926),  FliA  (GI:  1788231),  RpoS  (GI:  1789098),  RpoH 
(GI:  1789871),  Feci  (GI:  1790746),  andRpoN  (GI:  1789594) 
on  GenBank  database  (National  Center  for  Biotechnology 
Information  (NCBI),  USA)  were  used  as  queries,  and 
BLASTP  searches  [21]  were  performed  using  the  in  silico 
Molecular  Cloning  Genomics  Edition.  NCBI  conserved 
domain  (CD)  searches  [22]  were  performed  using  NCBI 
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Figure  2:  Circular  representation  of  the  Rothia  mucilaginosa  DY- 
18  genome.  Coordinate  markers  around  the  outside  of  the  circle 
are  in  base  pairs.  The  color-leveled  circle  indicates  putative  ORFs 
annotated  by  CRITICA  and  Glimmer  2:  ivy  green,  alternate  start 
codons;  light  violet,  energy  production  and  conversion;  blue,  cell 
division  and  chromosome  partitioning;  pink,  amino  acid  transport 
and  metabolism;  salmon,  nucleotide  transport  and  metabolism; 
light  pink,  carbohydrate  transport  and  metabolism;  bright  red, 
coenzyme  metabolism;  red,  lipid  metabolism;  teal,  translation, 
ribosomal  structure  and  biogenesis;  light  blue,  transcription; 
aqua,  DNA  replication,  recombination  and  repair;  dark  blue,  cell 
envelope  biogenesis,  outer  membrane;  brown,  membrane  protein; 
marine,  cell  motility  and  secretion;  bruise  blue,  post-translational 
modification,  protein  turnover  and  chaperones;  lilac,  inorganic  ion 
transport  and  metabolism;  orange,  secondary  metabolism  biosyn¬ 
thesis,  transport  and  catabolism;  light  orange,  general  function  pre¬ 
diction  only;  gold,  function  unknown;  grape,  signal  transduction 
mechanisms;  cobalt  blue,  transcription  regulator;  olive,  transport. 
The  inner  circles  correspond  to  a  GC  plot:  green,  below  average  GC 
content;  purple,  over  average  GC  content. 


server  (version  2.18).  Pfam  domain  searches  [23]  were 
performed  using  Pfam  program  (version  24.0). 

3.  Results  and  Discussion 

3.1.  Cell  Surface  Structures  of  DY-18.  Scanning  electron 
microscopy  revealed  that  DY-18  possesses  dense  meshwork¬ 
like  structures  around  cells  (Figure  1(a))  which  is  a  typical 
phenotype  among  bacteria  with  biofilms  [20,  24,  25]. 
ATCC25296,  a  type  strain  for  R.  mucilaginosa,  also  showed 
the  same  structures  around  cells  (Figure  1(b)).  In  both 
strains,  this  phenotype  was  maintained  after  repetitive 
passages  (data  not  shown).  Many  oral  bacteria,  such  as 
cariogenic  Streptococcus  mutans,  are  known  to  form  biofilms 
in  a  sucrose-dependent  manner  [26].  It  is  important  to 
note  that  DY-18  possesses  dense  meshwork-like  structures 
in  a  sucrose-independent  manner.  Considering  the  limited 
nutritional  condition  in  apical  periodontitis  lesions,  this 
behavior  might  be  important  for  this  organism  to  survive  in 
a  severe  condition. 


3.2.  Complete  Genome  Sequence  of  DY-18.  In  order  to 
understand  the  regulatory  system  behind  the  maintaining 
of  this  phenotype,  the  complete  genome  sequence  of  DY-18 
was  determined  using  a  combination  of  454  pyrosequencing 
and  paired-end  sequencing  [15].  A  total  of  193,024  reads 
from  standard  sequencing  and  83,610  reads  from  paired- 
end  sequencing  were  assembled.  Seven  scaffolds  containing 
64  contigs  were  obtained  with  99.71%  of  the  bases  having 
a  quality  score  of  40  and  above.  As  the  results  of  Gap 
closing  by  primer  walking  PCR  and  sequencing,  the  DY-18 
genome  sequence  contains  2,264,603  bp  and  had  an  average 
sequencing  depth  of  18-fold.  The  one  circular  chromo¬ 
some  had  a  high  GC  content  (59.6%).  No  plasmids  were 
found. 

The  sequence  was  annotated  using  the  CRITICA  and 
Glimmer  2.  The  entire  genome  contains  2,052  predictive 
protein-coding  genes  including  49  tRNA  genes  and  1 1  rRNA 
operons  (Figure  2).  Genome  information  for  the  chromo¬ 
some  of  DY-18  was  deposited  in  DDBJ/EMBL/GenBank 
under  the  accession  number  AP0 11540.  The  allocation  of 
main  gene  function  is  described  in  Table  1 . 

3.3.  Sigma  Factors  and  Two-Component  Systems  on  DY-18 
Chromosome.  As  described  above,  the  biofilm  formation 
is  tightly  linked  to  the  stress  responsive  system  in  many 
bacteria  [16-20].  So,  we  bioinformatically  determined  the 
stress  responsive  genes  on  the  chromosome  of  DY-18. 
This  organism  possessed  two  sigma  factors,  one  of  which 
is  a  primary  sigma  factor  that  recognizes  housekeeping 
promoters  and  the  other  is  an  ECF  sigma  factor  whose 
promoter-binding  activity  may  be  regulated  in  response 
to  physiological  and  environmental  cues  (Figure  3(a)). 
ATCC25296  has  one  primary  and  two  ECF  sigma  factors 
(accession  number:  ROTMU0001  _0778,  ROTMU000U0903 
and  ROTMU0001_0927).  Furthermore,  the  gene  encoding 
the  ECF  sigma  factor  is  arranged  in  an  operon-like  structure 
with  its  cognate  antisigma  factor  gene  (RMDY18_07080) 
(Figure  3(b)). 

Sigma  factors  are  involved  in  promoter  recognition  and 
transcription  initiation  [27]  and  plays  key  roles  on  bacterial 
stress  response  systems.  In  Staphylococcus  aureus,  the  sigma 
factor  influences  the  expression  of  several  virulence  factors 
and  is  associated  with  an  increased  ability  to  adhere,  invade 
and  persist  within  host  cells  [28].  Moreover,  the  activity  of 
sigma  factor  positively  influences  the  expression  of  several 
cell  surface  proteins  [29].  These  results  suggest  that  the 
sigma  factor  plays  an  important  role  in  the  establishment 
of  chronic  and  intractable  S.  aureus  infections  [28].  In 
Pseudomonas  aeruginosa,  ECF  sigma  factor  AlgU  responds 
to  environmental  stresses  with  the  production  of  alginate, 
a  component  of  its  meshwork-like  structures  [30].  The 
antisigma  factor  MucA  is  a  negative  regulator  of  alginate 
production  that  sequesters  AlgU  to  the  inner  membrane 
[31].  In  this  study,  we  attempted  to  find  genes  on  DY- 
18  genome  highly  homologous  to  well-known  biofilm- 
associated  genes  of  S.  aureus  and  P.  aeruginosa.  But  our  trial 
failed  to  hunt  out  homologous  genes  with  certain  accuracy. 
Instead,  BLAST  analysis  of  the  RMDY18_07080  gene  showed 
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(Primary  sigma  factor) 


7090 

(Additional  sigma  factor) 


Sigma  Antisigma 


7100  7090  7080  7070 


(b) 


Figure  3:  (a)  Domain  architectures  of  two  sigma  factors  on  the  genome  oiRothia  mucilaginosa  DY-18.  Existence  of  sigma-70  regions  1.2, 2,  3 
and  4  in  a  sigma  factor  7350  encoded  by  RMDY  18_07350  was  predicted  by  the  Pfam  search,  indicating  a  primary  type  sigma  factor.  Another 
sigma  factor  7090  encoded  by  RMDY  18_07090  was  predicted  as  an  ECF-type  sigma  factor  in  which  sigma-70  regions  2  and  4  were  detected. 
Proteins  are  represented  as  thick  horizontal  lines  and  their  domain  structures  are  shown  by  filled  rectangles.  Numbers  above  the  rectangles 
indicate  the  domain  names  (Sigma-70  regions  1.2,  2,  3,  and  4).  The  line  indicates  100  amino  acids,  and  all  the  structures  are  to  scale,  (b) 
Genetic  organization  of  the  RMDY  18_07090  (7090)  and  flanking  coding  sequences.  The  RMDY  18_07090  is  arranged  in  an  operon-like 
structure  with  its  cognate  antisigma  factor  gene,  RMDY  18_07080  (7080).  Arrows  indicate  the  coding  sequences  and  their  orientations. 
Two  neighboring  genes,  RMDY  18_07100  (7100:  lysophospholipase  LI)  and  RMDY18_07070  (7070:  5-enolpyruvylshikimate-3-phosphate 
synthase),  of  this  sigma-antisigma  region  were  not  regulators. 


that  the  antisigma  factor  of  DY-18  has  a  motif  respond¬ 
ing  to  oxidative  stress.  In  Mycobacterium  avium,  addition 
of  hydrogen  peroxide,  a  known  stimulus  of  the  oxida¬ 
tive  stress  response  resulted  in  elevated  biofilm  formation 
[32].  Alkyl  hydroperoxide  reductase  genes  (RMDY18_18480 
and  RMDY18_18490)  were  also  identified  in  DY-18.  As 
for  a  future  study,  we  will  examine  whether  these  oxida¬ 
tive  stress  response  genes  affect  the  biofilm  formation  of 
DY-18. 

Bacteria  also  exploit  the  two-component  systems  for 
controlling  gene  expression  in  stressful  environments  [33]. 
This  signal  transduction  system  consists  of  a  sensor  kinase 
together  with  a  response  regulator  and  is  highly  conserved 
among  bacteria  [33].  In  response  to  stimuli,  following  the 
autophosphorylation  of  the  sensor  kinase,  the  phosphate 
is  transferred  to  a  cognate  response  regulator  to  con¬ 
trol  the  transcription  of  target  genes.  In  Escherichia  coli 
strains  ATCC25404,  DH5a  and  MG1655,  the  QseBC  two- 
component  system  induced  by  autoinducer-2  is  known  to 
control  biofilm  growth  [34],  We  found  that  five  genes  in 
the  DY-18  chromosome  were  assigned  to  two-component 
systems  on  a  database  of  Kyoto  Encyclopedia  of  Genes  and 
Genomes  (KEGG)/PATHWAY  (http://www.kegg.jp/kegg/) 
[35].  Although  the  functions  of  downstream  genes  are  still 
obscure,  we  found  that  DY-18  possesses  orthologs  of  two- 
component  system  genes  responding  to  low  temperature 
(RMDY18_05700  and  RMDY18_05710),  cell  envelope  stress 
(RMDY18_06650  and  RMDY18_06660),  and  phosphate  lim¬ 
itation  (RMDY18-17310). 

Pathways  built  on  KEGG/PATHWAY  (http://www.kegg. 
jp/kegg-bin/show_organism?org=rmu)  also  suggested  that 
there  are  several  genes  relating  to  biosynthesis  of  sugars.  The 
pathway  included  the  genes  involved  in  the  synthesis  of  man¬ 
nose  (RMDY18_00300,  RMDY18JXB10,  RMDY18JHB20, 


RMDY18_00610,  and  RMDY18-17420),  which  is  a  major 
component  of  meshwork-like  structures  of  biofilm-forming 
bacteria,  such  as  Prevotella  intermedia/nigrescens  [20,  36] 
and  P.  aeruginosa  [30].  Furthermore,  twenty-eight  putative 
glycosyltransferases  were  found  on  DY-18  genome.  These 
genes  might  contribute  to  the  exopolysaccharide  biosynthesis 
resulting  in  the  formation  of  meshwork-like  structures  of  R. 
mucilaginosa  DY-18. 

The  genes  constructing  ABC-type  multidrug  transport 
system  were  also  found  on  DY-18  genome  (RMDY18_05690 
and  RMDY18_05680).  It  is  likely  that  the  biohlm  formation 
and  the  multidrug  transporter  contribute  to  the  therapy 
resistance  of  this  organism;  however,  further  thorough 
studies  on  both  systems  are  necessary  to  fully  explain  the 
mechanism  of  therapy  resistance  of  R.  mucilaginosa.  In 
immunocompromised  patients,  antibiotic  treatment  against 
R.  mucilaginosa  is  difficult  [37]. 

4.  Conclusive  Remarks 

As  far  as  we  retrieved,  this  is  the  first  report  of  the  complete 
genome  sequence  of  a  clinical  isolate  of  Rothia  from  an 
apical  periodontitis  lesion.  This  organism  has  a  unique 
meshwork-like  phenotype.  Additionally,  we  found  that  DY- 
18  possessed  only  two  sigma  factors  and  one  of  these  was 
the  ECF  sigma  factor  responding  to  environmental  stresses. 
Two-component  systems  also  exist  as  stress  responding 
systems  in  DY-18  genome.  These  systems  seem  likely  to 
be  crucial  for  the  biofilm  formation  of  this  organism  to 
survive  unfavorable  condition.  However,  further  studies 
using  microarray  constructed  with  the  data  obtained  in  this 
study  and  pathway  analysis  [38]  are  needed  to  figure  out 
how  R.  mucilaginosa  regulates  biofilm  structures  and  causes 
opportunistic  infections. 
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